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Abstract. Magnetohydrodynamic configurations with strong 
localized current concentrations and vortices play an impor- 
tant role for the dissipation of energy in space and astrophys- 
ical plasma. Within this work we investigate the relation 
between current sheets and vortex sheets in incompressible, 
stationary equilibria. For this approach it is helpful that the 
similar mathematical structure of magnetohydrostatics and 
stationary incompressible hydrodynamics allows us to trans- 
form static equilibria into stationary ones. The main con- 
trol function for such a transformation is the profile of the 
Alfven-Mach number Ma, which is always constant along 
magnetic field lines, but can change from one field line to 
another In the case of a global constant Ma, vortices and 
electric current concentrations are parallel. More interesting 
is the nonlinear case, where Ma varies perpendicular to the 
field lines. This is a typical situation at boundary layers like 
the magnetopause, heliopause, the solar wind flowing around 
helmet streamers and at the boundary of solar coronal holes. 
The corresponding current and vortex sheets show in some 
cases also an alignment, but not in every case. For special 
density distributions in 2D it is possible to have current but 
no vortex sheets. In 2D vortex sheets of field aligned-flows 
can also exist without strong current sheets, taking the limit 
of small Alfven Mach numbers into account. The current 
sheet can vanish if the Alfven Mach number is (almost) con- 
stant and the density gradient is large across some boundary 
layer. It should be emphasized that the used theory is not 
only valid for small Alfven Mach numbers Ma <C 1, but also 
for Ma ^ 1. Connection to other theoretical approaches and 
observations and physical effects in space plasmas are pre- 
sented. Differences in the various aspects of theoretical in- 
vestigations of current sheets and vortex sheets are given. 



1 Introduction 

Many structures in magnetospheres or in the solar corona 
are often described by quasi-magnetohydrostatic sequences 
of equilibria, e.g., |Romeou and Neukirch' (1999), Schindler 



and Birn ( 1982[), Becker et al.^ (2001) or Wiegelmann and 



Schindler ( 1995 1. Current sheets are important for storage 
and release, i.e. dissipation of energy, leading, e.g., to mag- 
netic reconnection, see |Neukirch| ( |2005bl l, or |Solanki et al.| 
( 2003 1 for observational aspects of coronal heating. The for- 
mation of thin current sheets, is often done within the frame 
of the quasi-static approach, where the time dependent pres- 
sure on the boundary is prescribed. Via the law of adiabatic 
change of the flux tubes in the frame of ideal MHD magneto- 
spheric convection, the flux transport takes place, and forms 
highly structured current sheets. 

Another approach of the formation of current sheets is a 
parametric method, performed, e.g. in |Nickeler et al. ( 2006| l 
and |Nickeler and Wiegelmann] ( |2010| l. Here the parameter 
is the Alfven Mach number instead of the control parame- 
ters like, e.g., in ,Forbes and Priest ( 1995[ l and ,Romeou and 
Neukirch ( 1999[ l. These parameters describe constraints or 
boundary conditions, like magnetic shear, pressure or cur- 
rents. 

Instabilities of current and vortex sheets, although with 



very high Alfven Mach number, have been analysed in Baty 



et al. (2003 i. In this article grid-adaptive simulations show 
that even in very weak magnetic field regimes (Ma — 30), 
the large-scale Kelvin-Helmholtz coalescence process can 
trigger tearing-type reconnection events previously identified 
in cospatial current-vortex sheets. 



The paper of Eyink and Aluie ( 2006) explains how for dis 
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continuous fields a part of the magnetic flux can slip through 
the plasma by an averaged turbulent statistic fluctuation of a 
subscale electromotive force (EMF), i.e. an electric field in 
the comoving systems of the ions. Another important aspect 
of this paper is that current sheets and vortex sheets both have 
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to exist and 'intersect in sets of large enough dimensions' to 
get 'ideal' magnetic reconnection, i.e. 'magnetic flux con- 
servation may be broken in ideal MHD by nonlinear effects' . 
Such tangential discontinuities are basically current-vortex 
sheets. 

From numerical incompressible time dependent adaptive 
MHD-simulations by Grauer and Marliani ( 1998 , 2000 ), it is 
shown that vortex-sheets align with current sheets. The ini- 
tial conditions in this investigation are given by the so called 
Orszag-Tang vortex. 

In the present paper we analyse and explain some of the 
physical implications of the non-canonical transformation 
technique (as developed and discussed by, e.g. |Gebhardt| 
[and Kiesshiigl |1992[ [Petrie and Neukirch[ |1999| ) to calculate 
steady-state ideal MHD flows with current sheets. We anal- 
yse the feedback-circuits on magnetohydrodynamic forces 
that arise, enabling a relaxed, incompressible steady-state 
plasma flow along a magnetohydrostatic background field. 
We analyse the connection between vortex and current sheets 
in the frame of this mathematical technique, but also for gen- 
eral field-aligned flows. 

In section [2] we discuss the basic stationary MHD- 
equations and introduce a general solution for the incom- 
pressible field aligned case, based on magnetohydrostatic 
equilibria. In particular, we investigate the forces induced by 
the flow. Section [3] contains special cases of initial potential 
and force-free equilibria. In section]?] we investigate the re- 
lation between electric currents and flow vortices in general 
and provide analytic example solutions in section ]5] Finally, 
we discuss the impact of this work on physical objects like 
planetary magnetospheres and the solar coronal plasma. 



2 Transformation from magnetohydrostatic equilibria 
to stationary MHD configurations. 

Under the assumption of a finite Alfven Mach number 
the equations of stationary incompressible and field-aligned 
MHD are given by 



V-v 
p{v-'V)v 



0, v-Vp^O, 

jxB~Vp, 

MaB 



V-B^O, 



(1) 
(2) 

(3) 

(4) 



where v is the velocity, p the mass density, j = —'V x B the 
current density vector, p the scalar plasma pressure, and Ma 
is the Alfven Mach number. The assumption that the mag- 
netic field and the flow are parallel leads to the conclusion 
that the electric field has to vanish everywhere identically to 
fulfill ideal Ohm's law. 

Within this work we solve the equations by an 
transformation-approach developed in IGebhardt and 



stationary one. We show that the general solutions of the sta- 
tionary, field-aligned incompressible ideal MHD equations 
lead to a very close connection between current and vortex 
sheets. We also show how the forces are compensated by the 
used exact solution method. 

In the following we transform the force terms to enlighten 
the connection to the magnetic and plasma force terms. Un- 
der the assumption of field aligned stationary incompressible 
plasma flow we get the transformation equations (see |Nick-| 
|eler and Wiegelmann| ]20T0l [Nickeler et aT ] ]2006l for de- 
tails). Keeping in mind that the mass continuity Eqs. (JTh, and 
(]4]l for every Alfven Mach number Ma and the density p 
satisfying B ■ 'VMa — and B ■ Vp = is a necessary con- 
dition, the general solution of the system Eqs. (JB-Q in 3D, 
but also in 2D, is given by 

Bs 
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where the subscript S defines the original magnetohydro- 
static fields. 

Thus one has a recipe to construct field-aligned, incom- 
pressible flows along magnetohydrostatic structures: the first 
step is to find a magnetohydrostatic (MHS) equilibrium, 
which we label with the subscript S for static {ps,Bs)- 



For the mathematical method introduced by Gebhardt and 



Kiessling'(1992) the static equilibria is written with the help 
of Euler potentials (a,/3): B — Va x V/3 or, if the equilib- 
rium has some sort of symmetry (e.g. in z-direction) with 
the flux function A: B = WA{x,y) x e^ + Bz{x,y)ez- The 
next step is to specify a Mach number profile M^(q;,/3), 
or Ma {A) in 2D, depending on the Euler potentials or flux 
function, respectively. One should be aware that the function 
p as a function of the two Euler-potential is independent of 
the choice of the Alfven Mach number. They are basically 
two independent functions. 

The analysis of the forces deals with all the force terms in 
the Euler-equation of ideal MHD, namely 

p{v-'W)v=jxB-Vp. (10) 

The identities Eqs. ^-^ lead with the Eqs. (JT]i-(]4]i to 



p{v-'V)v - 



Ml 



Mo(1-M2) 



iBs-V)Bs 



(11) 



Kiessling ( 1992[l, which transforms a static solution into a 



This means that the new inertial force, induced by the flow, 
depends mainly on the magnetic tension of the magnetohy- 
drostatic field, and of course on a function of the Alfven 
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Mach number, i.e. the inertial force is a transformed part of 
the magnetic tension force of the magnetohydrostatic field. 

We are now going to analyse the right hand side of the 
momentum or force equation Eq. ( [T0| , where we can identify 
the Lorentz force by 

jxB = fL, + fL2, (12) 

where 

is X Bs 



f LI — 



f 



I- Ml 
1 MaBI 



L2 



:VMa. 



(13) 



(14) 



Mo(l-Afir 

Here we recognize that the initial Lorentz-force is enhanced 
by a factor 1/(1 — M^) and changed by a term depending 
on the gradient of the Alfvenic Mach number and the mag- 
netic pressure. But the the term with gradient of Ma is also 
influenced strongly by the strength of the Alfven Mach num- 
ber itself. If Ma ^ 1 the term with the gradient of Ma is 
enhanced by the (1 — Af^)^-term, having of course even a 
stronger influence than the (1 — Af^)-term in the nominator 
of the modified original Lorentz-force f li. 

The second term on the right hand side of the momentum 



Eq. ( 10 1 is given by the pressure force Vps transformed via 
1 



Vp = 



-Vps 
Ml 
1-Ml 



MlBl 
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A*o(l-Ml) 

Thus after final addition of the new pressure gradient to the 
new Lorentz-force on the right hand side of the momentum 
equation Eq. ( [T0| the only remaining force-term is the net 
force, namely the modified magnetohydrostatic tension force 

-Ml/(l-Ml)(Bs-V)Bs^p{v-V)v. (17) 

Mo 

The remaining term 



Vps-jsxBs 



(18) 



vanishes, because this describes the original static equilib- 
rium. As it can be recognized by the inertial force (^tension 
force) in Eq. ( [T7] i, the expression can become extremely large 
in the limit of Ma — > 1. 

This implies that the modified pressure gradient consists of 
four terms: the first is the static pressure, the second a term 
completing the new Lorentz force to the old static one, the 
third term removes the gradient forces concerning the Mach 
number gradient term in Eq. ([12]), and the fourth term dis- 
plays the final inertia or net force term of the flow, induced 
by the mapping method. 



2.1 Comparison between quasi-static and stationary 
approach 

If the gradient of Ma is very large and also Ma has a value 
close to one, as explained in the last section and following the 



last paragraph of the Lorentz-force Eq. ( 12 1, hence the com- 
plete Lorentz-force can become extremely large. The situ- 
ation in the quasi-static approach is different: Although the 
current density grows extremely in the forming thin current 
sheet, no strong Lorentz-forces occur, due to the null sheet 
of the magnetic field, ignoring a small normal component to 
the sheet. 

In the case of the stationary states even for flows with 
small or moderate Alfven Mach numbers, the magnetic flux 
density and therefore the magnetic flux is not essentially en- 
hanced, compared to the initial MHS state, but the newly 
generated current density is able to increase in order of mag- 
nitude. Additionally, filamentary or 'fractal' structures com- 
pared to the old magnetohydrostatic current distribution are 
created, i.e. a double-structure evolves, see [Nickeler and| 
Wiegelmann (20101, and also see Figs, ([l] [2]). The pres- 



sure gradient is enhanced and has to be compensated by the 
Lorentz-force as the additional inertial force due to the flow 
is of second order in Ma, see Eq.([TT]i for small Ma ^ 1, and 
therefore negligible, but not negligible for large Ma- The 
only price one has to pay here is the unchanged shape of the 
fieldlines, in contrast to the quasi-static approach, and the 
formation of strong shear flows, i.e. small scale vortices, as 
we will show in the next sections. 

The additionally occuring filamentary fine structures of 
the steady-state current sheets, i.e. thin current sheets. 



are therefore similar to results found, e.g., in Wiegelmann 
and Schindler ( 1995 1. These authors found that only for 



non-similarity solutions of the quasi-static equations double- 
structure of the current density evolves. The evolution of a 
double-structure is necessary for the occurence of thin cur- 
rent sheets. 

However, the method of thin current sheet formation and 
structure, used in this paper, is a different approach than the 
quasi-static approach, e.g., see Schindler and Birn ( 1982|l 



or Wiegelmann and Schindler ( 1995 1, where the pressure as 
boundary condition is applied and via the law of adiabatic 
change is equivalent to magnetic flux transport and therefore 
flux enhancement. This is in contrast to the approach of this 
paper, where at least for small Alfven Mach numbers the flux 
is not essentially enhanced. 

One interesting aspect of the velocity fields, being ex- 



plained and calculate d in|Wiechen and Schindler] ( |1988[ l and 
Schindler and Birn (1987|, is that the velocity component 



parallel to the magnetic field in the vicinity of separatrices is 
getting very large. Hence also in the quasi-static regime, the 
flow parallel to the magnetic field is of huge importance for 
the plasma dynamics. 



Nickeler & Wiegelmann: Current sheets and vortex sheets 



3 Special cases 

To see how the transformation technique affects the initial 
MHS-equiHbrium we will now discuss and analyse some 
properties of specially chosen MHS-equilibria and flows, i.e. 
Alfven Mach numbers. 

3.1 Flows with global constant Alfven Mach number 

3.1.1 Potential field as initial static equilibrium 

The most simple cases of field-aligned incompressible MHD 
flows seem to be represented by potential fields as initial 
MHS-fields and the restriction to constant Alfven Mach num- 
bers. We now analyse some combinations of these possibil- 
ities. Discussing jg = and Af^ =:const it can clearly be 



recognized in the modified Lorentz force Eq. (12i that the 



magnetic field is similar to the old magnetohydrostatic field, 
i.e. the static field does not change its relative field strength 
across the field lines. The magnetic field stays potential, no 
currents form. The gradient forces vanishes, i.e. 'Vps — 0, 
there is no Lorentz-force, and only the constant static pres- 
sure PS changes to an inverted shaped thermal pressure p, 
reduced by a fraction of the magnetic pressure, depending on 
Ma- One needs large enough offsets in that case, to avoid 
negative thermal pressures p. The net force or inertial force 
from the flow can then be interpreted as the gradient force 
of the plasma pressure in the magnetohydrostatic case times 
Af|/(1 - M^), which for V x B5 = is identical to the 
magnetic tension force. Without gradients of Ma no addi- 
tional forces can be generated. 

3.1.2 Initial force-free fields 



tion equations are non-linear. This leads to additional non- 
linear effects as discussed in the following. 

3.2.1 Initial potential field 

The special case of j^ = and Ma T^const is simple, but very 
interesting for easy modelling: one can choose a potential, 
e.g. asymptotical homogeneous field and map it with mul- 
tiple current sheets, as done in Nickeler et al. ( 2006[ l. The 
effect on the forces is that the Bernoulli pressure force within 
the new pressure force, i.e. — Vp — —'V{ps — pv^ jT) must 
be compensated by an additional Lorentz-force. As the f i,\- 
term in Eq. ( 12 1 is zero, as js — 0, only the force-term f i^\ 



in Eq. ( 12 1 can do that job. The reason is that only the gra- 
dient of the Alfven Mach number can generate currents and 
then induce additional Lorentz-forces. The net force or iner- 
tial force originating from the flow can then be interpreted as 
the gradient force of the magnetohydrostatic pressure times 
M^/(l — M|). The fact that Ma is non-constant produces 
currents perpendicular to the magnetic field, repeating and 
keeping in mind that, 

VMa X B5 V X B5 



Moj^ Ma 



Ma 



VMa X Bs 



(1 



'M\)'l- 



(19) 



The equation reflects that in the this case the generated cur- 
rents are perpendicular to the magnetic field but also to the 
direction in which the flow gradient is oriented. 

3.2.2 Initial force-free fields 



The Lorentz-force equation Eq. ( 12 1 shows also clearly that 
instead of the initially force-free state a force is induced by 
the flow gradient term VA/^. All other aspects are similar 
to the potential field case. 

3.1.3 Non force-free initial magnetohydrostatic equilib- 
ria 

With this trivial case that jg 7^ and Ma =const, we get so- 
lutions very similar to the initial MHS equilibria. The mag- 
netic field and the current are mapped by a constant factor 
l/-\/l — A/^, and therefore the Lorentz force is enhanced by 
the constant factor 1/(1 — A/^). The current sheet structure 
is conserved, and therefore basically these solutions repre- 
sent similarity solutions, as all terms, including gradients of 
the Mach number, vanish. 



While the first term in the first line of Eq. ( 19 1 provides only 
perpendicular currents with respect to the magnetic field, if 
the Alfven Mach number has strong gradients, the effect of 
the second term is only to enhance the existing magnetohy- 
drostatic force-free current, being parallel to the magnetic 
field, writing 

VAf A X Bs V X Ss 



MoJ= Ma 



Ma 



{\-M\fl^ 
VMa X Bs 



{l-Mlf' 
aBs 



(20) 

{i-Mir^ {i-Mir^ 

where V x Bs — aBs- This implies that a former force-free 
field without flows, now develops components of the current, 
producing forces, pointing into the direction of "VMa, see 
modified Lorentz-force in Eq. ([T2]i. 



3.2 Flows with non-constant Alfven Mach number pro- 
file 



4 Connection between vortices and currents in field- 
aUgned, subalfvenic, stationary ideal MHD 



The generic case of plasma flows is that the Mach number 
changes from one field line to the other and the transforma- 



Even if the notion of 'turbulent' and 'ideal' reconnection. 



see, e.g. Eyink and Aluie (20061 or other references, are 
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based on discontinuous solutions, also including a specific 
angle between vortex sheet and current sheet, the solutions 
presented here, represents a powerful method to find the limit 
from continuous, namely well known MHS -solutions with 
smooth current sheets, to almost discontinuous solutions, i.e. 
sufficiently small subscales of concentrated and strong cur- 
rent sheets. It is obvious that for v\\ B the 'curls' of both 
vector fields, i.e. their vortices, are directly connected, of 
course not necessarily parallel[j 



MaB 



XB 



Vxv = VXxB + XVxB. (21) 



As V X B can be split into one perpendicular component j± 
with respect to the magnetic field and one parallel compo- 
nent, j II , we can rewrite the curl of the velocity 



V X ■» = VA X B H (ill +j±) 



This implies that 

B-{VxB) = 







(22) 



(23) 



is a necessary condition for a MHD configuration with a non- 
vanishing current but without vortices in 3D (or 2D). This is 
necessary as the term with the gradient of A, being perpendic- 
ular to the magnetic field, cannot compensate the term with 
j||. For a 2.5D configuration the parallel part of the current 
can only vanish identically if everywhere (V^)^ can be ex- 
pressed as a function of A. Thus, a physical situation with 
currents, but without vortices it is necessary that the parallel 
current vanishes. We will give an example, to be precise the 
sufficient criterion in 2D (Eq. ([25]l). 

With the additional assumptions that the flow is incom- 



pressible we can simplify the relations (21 1, getting the vor- 
ticity or vortex (strength) $7, 



n =vx 



/MaB 



MaVp xB V X {MaB) 



MaBs X Vp 



VMa X Bs 



MaV X Bs 



/JI^{l-Mlf' 



(24) 



For potential fields ( V x i?5 = 0) the vorticity only has com- 
ponents perpendicular to the magnetic field, if Ma is con- 
stant, then the vorticity is only non-zero, if the density is not 
completely constant. Strong vortex sheets only occur if the 
density gradient is strong enough. 

In the case of force-free fields there are nonvanishing com- 
ponents of the vorticity, even if Ma is constant, while the 
current remains field-aligned. The second term in the last 



equation also enables to generate a component of the vortic- 
ity perpendicular to the magnetic field in the case 'VMa y^ 0, 
but it is, of course, also perpendicular to the current. 

In the nonlinear case with small Ma =const, the vortex 
can be large, while the current is not extremely enhanced by 
the flow. Necessary for J7 = 0, e.g., in 2D with V x Bg = 
is, with apostrophes denoting derivatives with respect to the 
flux function A, 

VMa X Bg Ma VpxB 



V'(mop)(1-A/1) mo 2p3/2 
M'^{A){VAxBs) ^ p'{A)Ma{A){VAxBs) 
v/(Mop)(1"M2)3/^^ 2p3/2^^(l_Ml)^/^ 

M'^jA) ^ p'{A) 

Ma{A){1-M\{A)) p{A) 



P^Po 



Ml 
l-Ml 



1/2 



(25) 



The density p and the Alfven Mach number are usually inde- 
pendently conserved quantities along the field lines, i.e. both 
p — p{A) and AIa = Ma{A) can be freely choosen. One 
recognizes that this cannot be fulfilled, if the vorticity should 
vanish. Here, either the density is a specified function of 
the Alfven Mach number, i.e. p — p{Ma{A)), or vice versa 
Ma = Ma{p{A)). IfVxBs^G than a solution is only 
possible if again relation Eq. (J23]l holds, i.e. the current must 
be perpendicular to the magnetic field. Concerning the cur- 
rent equation, Eq. ( [T9| it can be recognized that this neces- 
sary condition Eq. p3| l must also be fulfilled, if one wants 
the steady-state-current to vanish. In the 2D case the current 
distribution and the magnetic field are then basically ID. This 
could, e.g., lead either to field lines being circles or straight 
field lines in 2D. 

In general (3D), if Ma is non-constant and Bs is a non- 
potential field, then the current can only vanish if Js|| = 0. 

If Ma is constant, or its gradients can be neglected, and 
the corresponding magnetohydrostatic field is potential it is 
possible to have no currents, or current sheet, but with the 
gradient term of the density it is possible to generate a vortex 
sheet in Eq.([24| at the same time. 



'For global constant MaI ^ they are strictly parallel. 



4.1 Sequence of equilibria towards almost magnetohy- 
drostatic configurations 

To analyse the connection between current and vortex sheets 
without loss of generality in 2D, we search for an ansatz for 
the Alfvenic Mach number sheet. We use an ansatz where 
the Alfven Mach number Ma is depending on one flux func- 
tion, but basically this view can be transferred to the prob- 
lem in 3D with two flux functions, i.e. Euler potentials. 
We present the principle how to construct a series of peaked 
function for the Alfven Mach number, sharply separating re- 
gions with flow from regions with almost no flow {Ma <C 1). 
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The Alfven Mach number should therefore be represented by 
a function series depending on a smallness parameter £ 7^ 0, 
which tends to zero, to mimic an approximately static equi- 
librium. This has relevance for a lot of situations in which the 
complete or large part of the domain have a basically quasi- 
magnetohydrostatic character, with some regions showing al- 
most no flows along the field line, e.g., Neukirch ( 2005 a[ l and 
[Romeou and Neukirch] ( |1999| l. The relevance for the case 
Ma ^ 1 is valid especially in the case of a low /3-plasma, as 



l>/3: 



P/P 



2 Ml 



Byi2^io) By{2^iop) 7 m| 



(26) 



7 being the polytropic exponent and Ms the usual Mach 
number, i.e. 



P 



- const , 



4: 






(27) 



Approximate incompressibility is plausible for Mach num- 
bers Ms <C 1, and hence the condition AIa ^ Ms ^ 1, to 
guarantee a low /3-plasma, is even more plausible. In fact. 
Ma ^ 1 is possible for regions in which the plasma /3 is much 
larger, e.g., in the vicinity of (degenerated) separatrices, be- 
ing null sheets of the magnetic field , or regions where the 
magnetic field in general is weaker. 

Taking a normalized flux function A and a scaling factor e 
into account, a very general ansatz for Ma can be written as 
the following sequence of functions 



Ma =±\Ma 



±l + tanh(f) 



lim M^A =0 

e-fO 



dMA{A) 
dA 



M'a{A) <oo 



^ M'a ^ e'^-^ cosir^ 
MaM', =e 



^A^'^^A 



.2n-l 



cosh ^ I — 



sinh(f) 
cosh^ff) 



(28) 
(29) 

(30) 
(31) 

(32) 



i.e., the terms with Ma and iW^ are bounded for e 7^ 0. A 
divergence of these terms and therefore the parametric for- 
mation of a vortex sheet can be recognized in the second row 
of the vortex equation, Eq. ( 24 1 if we consider lime -^ 0, and, 
of course, taking into account that VMAiA) = M'^{A)VA. 
For small Alfven Mach numbers, even if we do not only take 
potential fields, i.e. V x Bs — 0, into account, the domi- 
nant and important terms are that term in the vortex and cur- 
rent equations , that include at least the VA/^i-term. We can 
distinguish different regimes for the exponent n, which are 
listed in the following: 

- (i) For n > 1, the current sheet and the vortex sheet con- 



£ ^- 0. The reason is that the AIa M'j^ term converges 

with £2«-i=e2(i+|A-|)-i^^i+|5| torero. 

- (ii) For n^l the vortex sheet and the current contract, 
the amplitude of the vortex sheet stays finite, the cur- 
rent sheet amplitude converges to zero, also due to the 

M^ Af^-convergence, reading ex e. 

- (iii) In the range 1 > 77, > | the vortex sheet and the cur- 
rent contract, the vortex sheet amplitude diverges with 
£-(!-") for £ — > 0, the amplitude of the current sheet 
converges to zero with MaM'j^ c>c £1*1. 

- (iv) For n— ^ the vortex sheet and the current contract, 
the vortex sheet diverges with £^ 2 , the amplitude of the 
current sheet stays finite, converging with MaM'a cx £". 

- (v) For < n < i the vortex sheet and the current sheet 
contract, the vortex sheet (with cx £"^''^1 ,5 > |<5| > 
) and the current sheet amplitude diverge Af^M^ oc 
g-jl'^l; for n = the Alfven Mach number stays finite 
and therefore cannot converge to a quasi static case. We 
therefore do not take values of n < into consideration. 

These different cases imply that only for a special chosen 
funcion of the Alfven Mach number current and vortex sheets 
simultaneously appear, namely for the case of the exponent 
n < 1/2. For the other cases it is at least necessary that n < 1 
to generate a singular vortex sheet in the limit £ — > 0, but 
without getting a current singularity. Case (v) guarantees the 
existence of both, a current sheet and a vortex sheet. 

We did not discuss the influence of the V/o-term in the 
vortex-equation, Eq. ( 24 1, very intensively. The influence 



can also be very important, although it is with a cx A/^Af^- 
dependence of the same strength and quality as the depen- 
dence of the current density. 

4.2 General view on the connection between currents 
and vortices 

To clarify the structure of this connection between vortex 
sheets and current sheets from a different view, we introduce 
the abbreviations Ga = C^Ma x -Bs)/(1 - Af^)^/^ and 
Js ~ Js/{^ — M^y/^ and can rewrite the vortex-equation 
Eq. (|24]) and the current-equation Eq. ([19]) 



Ma 
1 



G, 



Js 

Ma 



n = ^GA + 
Vp Po 



Js 



MAp- 



/p 



n 



P " Ma 
2 p,o 



VpxB- 



V~P 



-G, 



-VpxB 

Ma 
Pop 



(33) 



(34) 



Ma 
( 

Mo 



(35) 



The terms of order O(Af^) can be neglected for small val- 
to zero, as can be seen by taking the limit for Eq.p2]i for ues of Ma- Then one can recognize that the most important 



tract with l/£ in Eq. ( 32 1, but their amplitudes converge 
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term, the zeroth order in Ma, is the G/i-term, followed by 
the gradient term of the density and the term with the current 
density as terms of first order in Ma- 

But even for the general case Ma ^ 1 the Gyi-term has 
the strongest influence on the generation of vorticity, as it 
scales only with p^^/^ in contrast to the j'-term with a 1/ p- 
dependence or the p^'^/^-dependence of the Vp-term. Only 
if the scale of the density gradient is much smaller than that 
of the Alfven Mach number the gradient term of the density 
can dominate. The G^-term is also favoured, as for values of 
Ma close to one the (1 — Af^)'^/'^-dependence enhances the 
influence of the VA/^-term on the vorticity stronger, than 
the (1 — il/^)^/^-dependence of all other terms in the vor- 



ticity Eq. (35 I. The nature of this transformation technique 
shows that these steady-state current-vortex structures are re- 
lated to static equilibria, and therefore the vortex structure 
and its scales governs the structure of the currents and cur- 
rent sheets. 



5 Examples 

As an example we investigate an initial 2D-equilibrium in the 
form 



A{x,v)=x^ + {jj2^l)y^ 



(36) 



which reduces to a potential field solution for a vanishing 
current j^ as shown in figure [Tlpanel A{x,y). We describe a 
Mach-number-profile in the form 



Ma{A). 



MAr 



1 — tanh- 



A 



(37) 



where Mai 



- 0.5 and d = 0.2 are free parameters, see top- 
left panels 1 showing Af^(A) and — -§a~ for the profile and 
its derivative, respectively. This introduces a plasma flow on 
some field lines and a moderate flow gradient perpendicu- 
lar to the x-point separatrices as shown in the top-right pan- 
els Ma and ^^ . As a consequence the initial current free 
equilibrium develops a current and vortex-sheet aligned with 
the separatrices (middle-right panels J^ and 51, respectively). 
The bottom-right panel ^f shows the ratio of these quantities. 
In figure |2] we use the same mathematical form, but 
provide a significantly steeper flow gradient d — 0.2/16 
along the separatrices , but a lower maximum flow velocity 
-^Amax = 0.25. The panels in figurel2]have the same mean- 
ing as in figure [T] We observe, that both the current sheet 
and vortex sheet increase in magnitude, but the vortex-sheet 
is stronger affected, as the ratio jf shows. 



6 Conclusions 

Numerical incompressible time dependent adaptive MHD- 
simulations by Grauer and Marliani ( 1998| [2000 ), show that 
vortex-sheets align with current sheets. Although it is not 



clear why this alignment occurs within the approach of these 
simulations, we could show that for steady-state, incom- 
pressible field-aligned flows this alignment appears naturally 
in a non-trivial way: even if the Alfven Mach number is 
non-constant, and non-trivial terms appear, resulting from the 
used non-canonical transformation mechanism, the resulting 
microscopic or thin structure and decrease of the amplitude 
of current and vortex strength show a clear alignment. The 
reason is that in both representations of current and vorticity 
the gradient of the the Alfven Mach number Ma generates 
the thin structure of vortex and current sheet. 

The larger the gradient, i.e. the smaller the lengthscale, the 
stronger the vorticity and the current. The difference is only 
that the current generating term is multiplied by Ma, i-e. in 
the case of approximative magnetohydrostatics Ma <C 1 the 
vortex is stronger enhanced than the current. As the scale for 
the flow, inducing the current sheet, must be very small for 
small Alfven Mach numbers, the solution can show strong 
vortex sheets but not so strong current sheets. 

It is also noteworthy that the function p can be chosen 
freely as a function of the Euler-potentials in the frame of this 
transformation technique. The free choice leads to strong in- 
fluence of the vorticity on the density gradient. This may be 
important for regions across magnetic boundaries, e.g. hel- 
met streamer regions, the heliopause or other astropauses, 
the magnetopause etc. where inner and outer field line 
and plasma regions are separated by separatrices or magne- 
topauses, separating regions of higher density from regions 
of lower density. Thus it is very difficult to get current sheets 
without vortex sheets, as the used non-canonical transforma- 
tion method automatically has to produce shear flows, to cre- 
ate filamentary strucured current sheets. These shear flows 
in fact imply vortex-sheets, hence generating current sheets. 

On the other hand for special configurations it is also 
possible to create a current sheet with a vanishing vortex 
strength, in 2D, where one necessary condition is automat- 
ically fulfilled, namely that the parallel current with respect 
to the magnetic field vanishes, i.e. jy = 0. But a special den- 
sity distribution must also be included to be sufficient: nor- 
mally the both free and independent integrals Ma and p are 
then necessarily dependent functions of each other in 2D, i.e. 
P^p{Ma{A)). 

The velocity fields calculated from the quasi-static ap- 



proach, see Schindler and Birn ( 1987 1 and Wiechen and 



Schindler ( 1988 1 show high bulk velocities parallel to the 
magnetic field, similar to the pure field-aligned flow of our 
presented calculations. There is also an analogy to the fine 
structures of current sheets in both approaches, namely the 
filamentary structure formation, i.e., the formation of a peak- 
ing, almost singular current, in time-dependent, quasi-static 
theory and the filamentary structure formation of the current 
sheet due to small scale gradients of Alfven Mach number 
of our stationary approach. The used technique to calcu- 
late exact solutions of ideal MHD from known magnetohy- 
drostatic solutions is very important to understand the struc- 
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Fig. 1. Transformation from an initial potential field with a smooth Mach-number profile (see equationlsTlwith A/Amax = 0.5 and d = 0.2.) 
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Fig. 2. Transformation from an initial potential field, similar as in figure [T] but with a much steeper Mach-number profile (see equation |37] 
with MAmax = 0.25 and d = 0.2/16.) 
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ture of current-vortex sheets also regarding eruptive space 
plasma processes like magnetic reconnection and following 
magnetospheric substorms, eruptive flares or coronal mass 
ejections. Even if anomalous viscosity is damped or mini- 



mized Ahmadi and Hirose ( 1981 1 and shear flows in specific 
situations lower the reconnection rate Karlicky ( 1989[ l it is 
important to generate smaller scale dissipation, see e.g., Baty 
|eraLl ( |2003] l. 

'Flare observations should be used to investigate whether 
the large-scale vortical flows, required to sustain the maxi- 
mum viscous dissipation rate, are indeed present in the flar- 
ing solar corona.', as theoretical investigations on steady re- 



connection and dissipation show, see, Craig and Litvinenko 
( |2009] l. 

Hence we expect that in the vicinity of current sheets close 
to separatrices, it is possible to observe such strongly sheared 
flows. From that point of view it would be interesting to 
search for small scale shear flows in the vicinity of thin cur- 
rent sheets to compare this with initial configurations in ac- 
tive or eruptive regions in the solar atmosphere or in magne- 
tospheres before magnetic (sub)storms. The involved scales 
may admittedly be extremely small. 

It is also noteworthy that in adaptive mesh refinement anal- 
yses of magnetic reconnection, the 'fractal' or better filamen- 
tary structure of current sheets can also be found and plays a 
very important role in fragmented reconnection, i.e. current- 



layer fragmentation, see e.g. Barta et al. (2011 1 and Barta 
letaLldSoTOl i. 

It would be of great interest to investigate steady-state 
current-vortex sheets for general field-aligned flows. 
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grant of the Max-Planck-Institute in Katlenburg-Lindau. 
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